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ABSTRACT 


The effort here was concerned with initiating a long 
term experimental study to develop the method of laser 
aii scattering as a plasma diagnostic technique. This 
was ane with a ruby laser and a grating spectometer cou- 
pled with sensitive detectors, at the steady state plasma 
facility at the Naval Postgraduate School. Although at- 
tempts to detect the extremely low levels of scattered 
light failed, excellent progress was made in determining 
existing apparatus limitations and in recognizing which 
system parameter may be adjusted in order to detect the 
scattered light; this will enable the future experimenter 


to progress with a minimum of overlap. 
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I. INTRODUCTION 


Phas thesis" isMintended as av*detailed report of the 
beginning of a long-term experimental effort to develop 
Thomson scattering of laser beams as a plasma diagnostic 
technique. The technique employs the particle-like behav- 
lor of photons in scattering off plasma electrons. In 
general the background radiation from an energetic plasma 
dominates making it impossible to distinguish the scatter- 
ed light. This problem may be somewhat alleviated by the 
use of an extremely intense, sharply tuned source in con- 
junction with narrow band detectors. The high power ruby 
laser is such a source and a grating spectrograph-photo 
multiplier tube arrangement is such a detector. Even with 
these improvements, a manifold of problems and complica- 
tions yet exist, not the least of which is careful planning 
of the geometry of the experiments and associated alignment 
of apparatus. Previous experiments have shown that laser 
scattering from electron beams of energy approximately 2 


kev and density approximately 5 x 10” em? yielded scatter- 


ing intensities of only about alo iat of the laser output. 
Though electron densities in the plasma under study are 


Significantly higher, this gives some idea as to the prob- 


lems presented in discerning the scattered light. 


It will be shown that the spectral distribution of the 
scattered light is in part determined by a dimensionless 
coherence parameter proportional to the ratio of incident 


light wavelength to the Debye length; 


— 

when @ is small, that is, when the wavelength of incident 
light is small compared with the Debye length the cross 
section is the appropriate Thomson cross section and the 
shape of the scattered spectrum is a Gaussian; the width 
is the Doppler width associated with electron motion and 
as such is a meaSure of electron temperature. 

When @ is large (@ >> 1) the spectrum has a cant iat 
line of the incident light frequency, Woe between two sat- 
ellite peaks at frequencies et oe where 25 is the plasma 
frequency. In this case the Doppler broadening of the cen- 
tral peak iS approximately proportional to the ion motion 
and gives information concerning ion temperatures. 

The research begun in the plasma laboratory of the 
Naval Postgraduate School has been under the cognizance of 
Professor F. Schwirzke and is financed under Naval Air Sys- 
tem Command Research Grant. The accomplishments made thus 
far and described herein will enable to the future experi- 
menter to progress with a minimum of overlap. The major 


obstacle yet to be overcome, and indeed the main problem 


in the development of this diagnostic method as a whole, is 
detecting the scattered light and distinguishing it from 
natural background radiation. Once this is done, the tech- 
nique may be rapidly and effectively developed. 

The diagnostic technique once developed, it may be 
employed to investigate turbulence ina laser produced 


plasma. 


II. METHOD AND THEORY 


A. ELECTROMAGNETIC THEORY OF LIGHT SCATTERING 
Following directly from the equation of motion the am- 
plitude of an electron oscillation caused by the electric 


field E' of the incident electromagnetic wave of frequency 


@ is; 
Oo 
_e E'(t) _ ed7E'(t) 
SC) a 2 = 2 (1) 
m Ww 4? mec 


Multiplying the displacement by the electron charge e 
the induced dipole moment is obtained; 
Zid 
- E 
| Se eles es) = ian — a — A = = €E' (2) 


An-me 


For the laser light scattering to be dealt with here, 
is a plane monochromatic wave cut off exponentially at 


time T = 1/y, the pulsetime constant. 


That is, 
E'(r,t) = E' explilk'-r + (w+ iy)t]} (3) 
p(r,t) = ecE'explilk'-r ee ce isp”) edhe (4) 


It is this time dependent Dipole moment (4) which produces 


the scattered radiation. 


10 


The. ehbectric field ey ‘Zat) scattered by the a electron 


at position rs @etected at k,.from thesHertzZ.fommdae.for 


dipole radiation, is 


E. (R,t) = £.(R.,t) expli(k'-r. + k'R. 5 
2B. (R, t) a : ) expli(k #s 3) (5) 
R. = R. - r. 6 
J oa re) O 
where R. is the distance between scattering center and de- 


tector 


ee) is the amplitude of the scattered field and is 





2 
(i nee Ge leep iia — yield (7) 
po R.2x2 J J = O 
J 
a 2 
Riess, Mel. © 22" 2) 
4? mec 
= 1 A A ‘ : 
POER.,t)e= = =.) — (4. X(a.X E Jex 1W t-yvt 8 
£,(R.,t) ca la aaa aan 
Re 
where Sa is the unit vector in the Ra direction 
5 
i is measured from a reference point, 0, in the scattering 
volume. 


If the dimensions of the scattering volume are 


small compared with R, then Ro may be replaced by Je and 


Li 


£ = f for all electrons in the plasma volume. When deal- 
ing with the phase, however, more care must be taken in 


approximations. The following substitution may be made for 


Y see ie 
ry ze 
———_———— _~(Q'r- 
R. = /(R-r.) 2 = J (R2-2R-r tr .7 = eee to soee (9 
~J a ~ 7) eae 
Writing the scattered wave vector k" as having approximately 


the magnitude of the incident wave vector k' and direction 


R, that is, 


i 
k =k! Pay (10) 


E.(R,t) = £(®,t,rA8)exp(ik'R + ik-r.) (11) 


where k is the scattering vector. That is 


" ; he ; 8 ‘is 
k =k -k , k = (5 )sinS, k! wk = 4 (12) 


The scattering angle 8 is between k and a The direction 
is of course determined by the incident laser beam and by 
the position of the detector. 


In 90° scattering, 
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Furthermore }'~ X.=. The wavelength difference between 
A4' and A" is usually small and may be neglected, such that 
the symmetric of the scattered spectrum is undisturbed; if 
w is the frequency of the incident light and W is the fre- 
quency shift of the scattered light, T(W +W) a T(W -w) "no 
drift velocities exist. 

The phase of the scattered light depends on the posi- 
tion of the scattering center by the factor iat es «, ties 
becomes important when the fields scattered by a number of 


centers are superimposed at R. Summing, over all electrons 


in the scattering center, one obtains 


N N. 
; -ik'n Vik'TS 
B"(R,t) = )E.(@,t) = £(8,t)e ye (13) 
an Tal 
The scattered intensity is equal to the time averaged 
41 2 
Poynting vector we | a Introducing the incident light 


intensity as 


2 
tate (14) 


oO ANT 


one obtains 


r(@.4) = 10) = 018) <1), Thee s (15) 


EG 


Where 0@(8), Thomson's scattering factor, is equal in essence 
to Local. Introducinguewand.c. jethe polarization vectors 


of the incident and scattered light, 


E' = e'| B'| , E = e"| 5" | : Ree" _—" (16) 
2 2 
f i 2 
Gung) = Fal = = Ty e - (€e') rs (17) 
R 





for linearly polarized light; and for natural incidents 


ight; 
Pee 2 
“o(8) = 1, (£5) {1 + cos’ 9) (18) 
mc 2R 


where the 1 in parentheses represents the contribution 
from natural incident light and cos76 represents the con- 
tributed polarized parallel to the plane of scattering. It 


is easily seen how (17) complies with this. 


1B SPECTRAL INTENSITY DISTRIBUTION OF SCATTERED LIGHT 

For the most part, it is the electron density fluctu- 
ations in the plasma which are responsible for the scatter- 
ing of light; that is, according to Huygen's principle no 
scattering can occur in a strictly homogeneous medium. The 
electron density n(r,t) is a funetionsef pesitiongand tiger 
and may be considered as a continuously varying quantity or 
aS a superposition of delta functions as circumstances 


warrant. 
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Imagine the plasma confined to a large cubical volume 
V= _° and write the electron density as a spacial fourier 
series. 


igexr 


n(r,t) = ni(t)e = (19) 


<q ih 
ha hq 


Where q is the wave vector satisfying the cyclic boundary 


conditions 
27h} 
q; = ah h = 0,41,#2 , h. = X,Y,2 (20) 


The fourier coefficient ~_— be considered as the ampli- 
tude of the longitudinal plasma wave with wave vectors gq 
and may be written 


Ny - J mee) exp (= 1eie2a) Gr (21) 


for the microscopic representation and 


N 


ny = ) See | (22) 
jet 


for the microscopic representation. 
Usually, longitudinal oscillations in plasmas have com- 


plicated time dependencies and dispersion relations; further- 


more a number of frequencies may be present. The probability 


1 
As done in Light Scattering and = Instability in 
a Plasma by Otto Theimer. ltc, 
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of a certain frequency occurring is related to the mean 


Square of the amplitude ng () of a wave with wave vector 


q and frequency W. A long lived, relatively undamped 


oscillation will have a sharp maximum in Ng (we) . A con- 


tinuous function of the frequency Ng () is the Fourier 


transform of Ng (t) with respect to time and may be written 


iAP Treeeral sora 


CO 


Dee = | ng 66) exp (-1Wt) dt (23) 


But the actual density fluctuation is not by itself of pri- 
mary concern in scattering theory; rather, the amplitude 
of the polarization waves Pak) induced by the incident ra- 


diation is of interest. [Not to be confused with the dipole 


moments of equation (2) ] 


Bt) = € nee (t) (2a) 


Where € is given by (7a) and is the polarization of the 


free electrons acalm, ae, is the electric sre ld Omeseme 


incident laser pulse and has the time dependence 


iy (Ey, 


E exp[ (iw -y) tJ ’ t > 0 


(25) 
E (t) 


H 
S 
~ 
(Tt 


< 0 


That is, the laser is fired at t = 0 and is essentially 
overafter T == - Thus it is seen that the polarization 


waves have a frequency dependence (+0) the amplitude of 
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which decreases exponentially with the laser pulse, as 
exp(-yt). That scattered light will of course have the 
same time dependence, and the scattering waves may be iden- 
tified as a frequency shift in the scattered spectrum. 

Thus the finite laser pulse may be attributed (at 
least mathematically) to the plasma waves and A (w) may be 


written as 


fa 9) 


ng (w) = k ng(t) exp [-(iw + y) tlat (26) 


Considering the total scattered intensity, the time aver- 


age value of the function nm, (t) is taken 


Ppt) * exp (-2yt) dt 
‘e) 


I 


<|ny.(t) |“> 


| exp (-2yt) dt 


- ( 2) 


2Y | In, ce) |? exp (-2yt)dt 


— © 
, 2 
By Parseval's Theorem from (26) and (27) it follows 
that 


f=) 


C[n,,(t) |“) = x | _ yf) |* aw. (28) 


Clearly, additively combining the square of the amplitudes 


y a. | 2 
ie) J ny (w) | gives the square of the amplitude C] ny, (e)| par 


2 alin: 
In Bronstein and La. Semgndjajew, Taschenbuch d. Mate- 
matik,Verlag Harri Deutsch, Frankfurt/M.,1965,va5 pp.404-474., 


ly 


Thus each wave of wave vector k and frequency | contributes 


to the total scattered intensity by an amount given by 
= y ) - Cees 
I(w ,k)aw = 0(8) CF ) | n,(w) [aw | “25 (29) 


And as already stated, each frequency shows up as a fre- 
quency shift in the scattered spectrum. 

So it has been established that each electron becomes 
Hertz oscillator, when induced to oscillate by the inci- 
dent light, and itselfemits the light which is the scatter- 
ed radiation. If all electrons oscillate in phase with the 
incident light, the emitted (scattered) radiation will be 
coherent in time. If, however, the electrons have complete- 
ly arbitrary spacial distribution, coherence is destroyed. 
But electrostatic interaction prevents the electron distri- 
bution from ever being completely arbitrary, and it is from 
partial spacial coherence that so-called collective scat- 
tered radiation is observed. The degree or intensity of 
coherence is dependent on the ratio @ between the wave- 
length and the characteristic Debye shielding length for 


the plasma. That is, 


i 
os r Z KTs = = 
(30) 


Where X is the wavelength of the incident light, Xp is the 
Debye length and 9 the scattering angle. For @ << 1 scat- 


tering is incoherent; for @ >> 1 spaeaal coherence 
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determines the shape of the scattered spectrum. Perhaps a 
more physical and qualitative (if not mathematically rigor- 
ous) understanding may be had of this concept if one imag- 
ines the Debye spheres as integral units. Imagine the elec- 
trons as tied to the ions they shield. Within the Debye 
sphere, the electron distribution is arbitrary; that is, 
incident of wavelength less than the radius of the Debye 
sphere will oscillate electrons within that sphere in dif- 
ferent phases, whereas incident light of wavelength greater 
than the Debye sphere will oscillate the Debye spheres, and 
thus all electrons within them in unison, in coherence. 

The incoherent scattering spectrum is a direct illus- 
tration of the electron velocity distribution, as will now 
be shown. As the scattering under consideration is inco- 
herent, the effect of a single electron may be looked at. 
It has been seen that the electric field of light scattered 
by the sa electron, by equations (7) and (11) is propor- 


tional to the phase factor expliw t = ik-x,J. Approximat- 


ing eae by an expression valid for straight line electron 
paths, 
Percy = ©. (Opa Ve ou 
i aA) ae We ’ (31) 


the scattered field has the frequency 


Ww = W + kev. v2) 
2) aaa aa 


and the scattered light is Doppler shifted by E-y, = ie 


(k = eae. The shift is two fold; that is the electron 


oS) 


is a moving source to the scattered light. Thus the total 
shift is kev = k'-v - k'+«v. F(v) now, is the normalized 
velocity distribution function with velocity v parallel to 
k. Corresponding to velocity v, the frequency is W = kv; 
the velocity interval dv correspond to frequency interval 


Gdwe-—sker, Thus, hemnunbereor ~llectePons With velocity, 


1S proportional to the scattered intensity with shift kv. 


That is, 
k 
F(v) = SE (33) 


where v is the velocity component parallel to k, I(kv,k) is 
the scattered intensity integrated over all frequencies. 
Again, Chis is Onivenalicuror Ge —— 1 

The topic now turns to the more complicated theory of 
longitudinal plasma waves. The attempt here is to calculate 
the mean thermal value of the square of the space time de- 
pendent Fourier coefficients | ny, (w) |?. The basic reference 
for this calculation is the sane by Salpeter.> The ampli- 
tude squares | n, (w) | * wild be solutions of the collisionless 
Boltzmam eae en: i.e. the Vlassov equation (See Appendix A 
for derivation of the Boltzmann equation). The collision- 


less Boltzmann equation, 


(34a) 


Shed 
ct |Fh 
3 |e 
Fh 
i 

(= 


+ v°V £ + 
— 8 


S 
Pte Godeeter, Dip Smeghec ) Qaeedo 28 (lOO Oe 
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is essentially a continuity equation for representative 
points in electron phase space. The increase in time in 
the number of electron points in phase space is caused by 
three factors: (1) The drift in points in position coordi- 
nates due to inherent velocity: (2) the drift in points 
in velocity coordinates due to electric and magnetic forces, 
and; (3) the spontaneous velocity change without position 
coordinate due to collisions. 

In the purpose of calculating | n, (w) |? the spacial 
Fourier coefficient see of the distribution function 
P(e, Vv, t) 1s Lee (q is the general wave vector). 


Corresponding to equation (22), 


(vdv) 

ove = ) exp[-ig-r.(t) J ; (34) 
J 

the sum being taken over those particles with velocities 
between vtdv. As dv ~ 0, the right hand side has one term 
im a particle’ of velocity v®exrsts and no't€rm TF"such "a 
Merticle doesnot exist. “fhus (34) represents a fine grain- 
ed “is'tribution furetrvon "so described because the density 
has but two values: zero in a vacuum and infinite at the 
position where a particle is represented by a mathematical 
point. Clearly, then, an integral over the velocity dis- 
tribution becomes identical with a sum over all the parti- 


cles: 
N 


| if vy bt) dy = ), expliger, (t)] (35) 
be 


2 


The number of terms in (34) is equal to the number of par- 
ticles with velocities in v and v + dv. If the physical 
properties that are influenced only slightly by the micro- 
scopic details of the distribution function are calculated 
and if the course grained (continuous) distribution is in- 
dependent of spacial position, the number of terms in the 
sum is the most important aspect of (34), and may be approx- 


imated by 





a = n(z, bev )o= — ), 9 (t) exp (ig) (36) 
q 


where F(v) is the normalized coarse grained (continuous) 
velocity distribution and n ft) is the Fourier coefficient 
of the particle density FG by (21) and (22). As the no- 
tation in (36) implies, F(v) is independent of time. This 
of course, is valid only if AF/F is small during the time 
1/Aw where Awis the practical spectral resolution. 

Without regard to the coarse or fine-grained structure 
oir Ea ylt)s the distribution function may be represented for 
ean electrons aS a Fourier series: 


Baraat ; _ 
Po(fevet) = > »s a Gyo (cis , Ss i,e (37) 
q 


Similarly, for the three terms in the Vlassov equation, 


he 


Of cate) 

Sit il sc v 

ae TV ) 8t exp (1q-X) (38) 
q 





Ze 


and for the second term 


S LV 
ag alls i ), BE gqy (t) exp (1q°r) (39) 


In the third term, only the electrostatic force e E(x, t) is 
considered where e_ is the charge of the particle at r-E(z,t) 
is written as a Fourier series, 


E(x,t) = (t) exp (igq-r) (40) 


E 
“a 


<i 
LQ [7] 


For a equation (36) is used, where F(v) is the only ve- 


locity dependent quantity. The third term then becomes, 











7. Ct e. y OF. (v) 
Be age = Soy DLE explig’ +z) Jn, vex (ig"-x) 1-—S— (4) 
S a mVq- Gi > 
S ct at 
Combining terms in the double sum for which 
qq? q"™=q ) g"=q- q' (42) 
_ of . “s Pp ehe (w) 
= Et Pa EPRELS 0) (REG 036) aurea (43) 
m_ ov mv" o ie s(q-q') = 3v 
OF (v) 
ees i 
= 5 Pz. exp (iq-r) = SV 
mV oq = 
S =, 
where 
i 
sq °s0.Aq' “s (q-a') 2) 
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The "Random Phase" approximation is now introduced. The 
order of magnitude of ae follows from (22). Remember that 


Lim n, =N >> ae - /N (45) 
Nwe 


be is the sum of N positive and negative term with values 
arbitrarily distributed from zero to one. The probability 
that the real and imaginary parts of oe have the values 


and ” respectively is 
ae 5 p y 


(Eg + Ny”) 
1 Sle Sl 


and the mean amplitude square value is 
<a, |=] [ agian (é.? + 1 2ywié_.n) = 2N (47) 
q ree emanate) q q’ g 


(4) and (47) explain the importance of (45). ES that is, 
q = 0 has no dominating position for the electric field 
because it has a mean value in a neutral plasma of zero. 
Thus (44) has one dominating term only. 


Writing explicitly, 


a 15; = ee [EN +f DEF stag") | (48) 
ata 
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The “random phase" approximation assumes that the sum can 
be neglected in comparison to the term EN." Adopting this 


approximation, equation (33) finally becomes 





a of sav fgts oF =) 
V Jomo (igr2)| pa ¢ i Me Ea ait m eq ov |= . 


(49) 


The Fourier representation is unambiguous, and thus (49) is 
valid not only for the sum of all general wave numbers gq 
erie also for severy single q. Writing (49) explicitly fer 
electrons and multiplying by the*electron charge e, the 


equation for the Fourier coefficient o. ah corresponding 
eUrv 


to charge density Foy? 


30 Z 
a en OF (v) 
ae + 1vegq 9 oo cae a = 0 (50) 





An analytical expression for E may be obtained from 


Poisson's equation 
° — = Ls : 
V -E ATE — AAS Oe) (51) 


Where 0 is the charge density created by electrons and ions. 


As previously done, let 


C{r;e) = ) a, (t) exp (ig-r) | (52) 
q — 


4 
Though the infinite number of terms in the sum makes 
it by no means really clear that the approximation is al- 
lowable. 
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Bacar 5 2d al SIM. alte ( 4.00) 


q 
) ig-E = 4 do ,» E_ = -4%10 =5 (53) 
q g q a belo; 
and q = k (as in equation (12)) 
SF oky ann.* OF, (v) 
— eee ad ir=) i, = 1 ( +a)  ——— (64) 


Whemey V_ and, E aregsthescomponents ofev and, Eepara Wlelatosls 


(k is along the z axis). The Fourier time transform of 


Coy is needed in the calculation of the spectral distribu- 
tion; this is written as in equation (20) 
foe) 
apy (H) = | Saycy(thexpl-(iw+y) tat (55) 


Integrating over all velocities, 


os) CO 


| Lore Qo. (4) = | 7. (t)expl-(iwty) elat (56) 


iI 


Dory (“ Satisfies partial differential equation obtained .by 
equating the transform on the right and left of (54). The 


following integral is needed 





role) 
ay expLl-(iw+ty)tl]dt = ~F yey (0) + CAW+Y) Gay (Ww) (57) 
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This is obtained by integrating by parts and using (55). 
Notice that the contribution from the upper limit oe, aw 
disappears due to exp(-yt). Using (57) in (54) gives 


4Tn 


OF 
Dany (#4) (Lety+ikv) = O oye (9) + 1 ( = 


e ) one 
Q.(a) ~—Z— (58) 





Semualtiplying by 1 





4m. ae (v) 
4 (0) -C = oo, te a 
ee - — me 
where 
(4 = (4) + 2, (4) (60) 


The charge density for the ions is obtained similarly, that 


is, 


OF. (v) 
4 (0)= at as 


male) - -@ - a +1LY (61) 





and the total charge density 


4m. oF (v) oF (v) 
Soy (0)~ (oon a, (uf Se = 


Toy 





, ~ dex 7 OFT (w) = -W - kv +1W 


(62) 
Integrating (62) over all velocities gives an expression for 


Q. (4), which may then be introduced into (59) and that 
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result integrated over all v. The final result would be 

the amplitude necessary for scattering Qa, (4) = - en, (4). 
Since the integrals may not be actually evaluated without 
choosing special velocity distributions, the following 


function is antweduced 





; (iy ) i" oF. (v) dv 


G_(w) = SS - 
S M 4 J 4 SV  —-Wkvtiy 


s = 1,e. 


(63) 
Furthermore, corresponding to (35) the velocity integral may 
be calculated as a sum if id, (0) is chosen as a fine-grained 
fae wont, - 


Then, by integration of (62) 


exp (-1kZ.. ) exp (-1kz.. ) 
ee 
-wrkv thy “W-k tly 
() =—> 7 SS Ee 
o).4¢) ie ) 1 “cea Ge (64) 
7 e 1 


Where Be AC WC Ve5 are position and velocity components 
of ions and electrons parallel to k at t = 0. 

Equation (64) may be considered the basic equation of 
the theory of longitudinal plasma waves. From this equation 
follow the formulas for scattering, dispersion and drift 
instabilities. 

The next step is to employ the theory of density 
fluctuations to light scattering ina Plasma. This requires 


the amplitude squares | Q, (4) jos and 1, (w) io these 


represent a time average and a time integral over all micro- 
yp 


scopic states assumed by the plasma between t = 0 and t = y* 
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Substituting (64) into (59) and (61) and integrating 
over all velocities using the same method by which Q,. (4) was 


obtained, lo_, (wl * may be calculated. 


(1-G. (W)) 2 exp (-ikz.) Y exp (-ikz.) 
ci. . sume 2 1 —— Cs | aoe 
ex (“= Te (w) -G, (@) | ~wkvity 8°) Titkvsiy | 
(65) 
and 
i 1-G_(&))2 exp-ikz.) u exp tikz.) 
0. (W= a_ts | "Se ae ead 
ek 1-G, (4) SC ? —-Ww-kv4+1y G, (4) ae 
(66) 
The amplitude squares lo, (wl? contain double sums of the 
ewe 
Z exp (ikz. Lo explik(z,-z.) J 
: ~wkv .+1Y =i? (—W-kv_.+1y) (-@kv_ ,+iy) (67) 
ej ej el 


as well as double products. Ion velocities and positions 
are indexed by £, whereas electron velocities and positions 
are indexed by j. Diagonal terms of the left (the single 
sum) are treated the same as off-diagonal terms. In the 
double sum, however, the diagonal terms behave differently; 
that is, the space dependent phase factor vanishes in these 
terms allowing the replacement of the sum with an integral 
over the velocity distribution: 


N 


) explik(z.-z,) ] “ F(v) dv 
reo - "| Syne (68) 


oO 


Since the denominator of the right hand integral depends 


only on Vee the integral is over the one-dimensional 


Zo 


distribution Fi(v), gotten by integrating Elite) over v,, and 
aes Rigorous evaluation of the integral is at best tedious. 
Beaicly, ~~ewever , be approximated more readily. The main 


contribution to the integral comes at resonance v = - — 
At resonance, the value of the denominator varies rapidly 
while the numerator varies slowly, permitting it to come 


out of the integral. 


— W 
en 
2 Se ee 2 2 vk 


, F (v) dv 
-° (Utkv) “+ y 


lim a 
yrO }=—-@ (W+kv) + ¥ 


(69) 


yielding an expression for the diagonal sum dependent on 
the time of observation — This, of course, is what is 
expected, since |, (a) ie determines the scattered radia- 


tion during T = 


= e 
6 

The off-diagonal terms behave somewhat differently. It will 
be shown that this sum may be neglected compared with the 
diagonal terms. The sum 1S independent of Y and may be 
neglected in comparison to the double sum (which depends 

on +) for small y. Consider the double sum over a factor 
explik(z,-z,) 1] and assume the denominator is replaced by an 
appropriate mean value. The result 1S a sum over N(N-1) =N7 


POsrttive and=negative terms. From the®same “random phase, 


the sum is of the order of magnitude N. Hence, 


explik(z, =Z.) f [Peleg Few avjavy 

ae ce eee) Os Se ae 
) ) (-a-kv, +i y) (-aekv. -17) S Jp COOK, FEY) CO “kv y-1Y) 
i) Lee J ce 
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The single integrals over ane and av, have linear resonances, 
which by symmetry are integrable; the numerator remains un- 
changed as the denominator changes sign. Thus the resonance 
contribution to the integral is zero. Therefore the single 
integrals are independent of y and may be neglected in com- 


parison to the diagonal terms proportional code. - “Onémthen 


obtains, 
2 2 —w yi2n cw 
eae Ce | 1 (2 mn LS (4) ee Ga (4)| Fs (ae) 
2 eee! = ae 2 
e | 1-G. (uw) -G, (a 
(71) 
y| 2 2p (fw yep ce 
LOCO my | 1-G (I FL G+ 16 (TFL GO nf 
zZ yk Z 
e | 1-6, (4) -G, (4)| 
Hay omy FOG +E 
e YE | 1-c.(w) -c, («| ? — 


Evaluation of these amplitude squares of course depends on 
calculating Gy and G. for special electron and ion velocity 
distributions. 


Of interest is the shifted Maxwellian distribution; 





: aie. ( wart) 7 ae | - are Ov w_)? | (74) 
3 S 
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where W. 1s the drift vele@en tyeot wpartme@les»,of eye S;, 


(S = 1,¢é) 
Clearly, 

OF (v) m 

Se Ley 8) ] Fe (733 
and 


(vtW_) FP. (v) dv 


2 
— ‘e'd - 
oo ae | eer ek vty) 


(76) 


where a. is given by (30) and G. (4) by (63), and where 
Bev) is a one-dimensional velocity distribution in a and 
i is the z component of We 


Introducing the Doppler shifted frequency, one has: 


W = W- kW (77) 
SW S 


Quantities of the unshifted (no drift velocity) distribu- 
tion will now be superscripted with a zero. 


Then, 


Loe) 





O 
3 (VWs ) Be 2 ican A 
s | Ss _o [-wW_ -k(v+w_)tiy] Cale.¥) (78) 
Sw S 
and 
” BL (v) dv ~ Bo (vaw_) d(v4w_) -W 
S - | s S ee ee 0 (. Sw 
© (utkv) 7+y7 -o fw +k(v+w ieeo eS SE 
Sw S 

(79) 


a2 


Using the Doppler shifted frequencies 


2u 

ws? = =e a s,s' = i,e (80) 
s,s 3 s,s a2 
Ss 


where ve = is the mean thermal velocity square of the par- 
e 


ticles of mass ™m and temperature To. and introducing the 


dimensionless frequencies 


ew ot: 
ee re 


it follows that (using (74)), 














| Te a a ae eee (y' ‘IF exp (- seat ae (i |? ae 
ny, (4) 16,0) 2 ey TE 
(82) 
and, 
cme L W 
yaeG@-@ ft. w -B-O =| © 
xl Fe a} ii ee 
where 
m hs mM, L. 
— a) GD (84) 


Furthermore, if the substitution 


u = Wey (85) 
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is made, and the integrand of (79) is expanded with the 


factor (-W-ku-1y), then, 


= O 
Cen = Ko | uF (u) (-@Ku-1Y) du 


a? . ¢e Gi) Secs as 
The imaginary part 1S an integral of the form of (69). 
Therefore, 
lim I(x) = 7% ox exp(-x*) (87) 
a0) 
Assuming yY small the real part of G(x) may be expanded in a 
power series, which, when evaluated, is equal to zero in the 
limit. Using the same arguments as with (70), the real part 
goes through zero at resonance and changes sign. Again, 
integrated over symmetric limits, the integral goes to zero. 
From toot let <| ul> be. the mean value of | ul; expand- 


ing the integrand of (86) for w>k<lul> in the form 


(o= 





Lin yes = age a a A” (32) 
y7O (-&ku) +¥y FG 
the following is obtained: 
a 3 eee 
lim Re(x) = —— (Ll 4+)3 23ers) (89) 
2 Z 4 
y70 
For w<k <Jul>, the substitution, 
=k we, v= - oO, au = - 2 (90) 
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is made and E> (u) 1s expanded. That is, 


“7 4q\¥ >, m 
PaO HE) = 6 Deiat Hey 2) 


v=0 
(99) 
Using this to obtain 
lim Re(x) = - a*[1-2x7 (1 - £ x? + = eee -..)]= 
veO 
2 
- a [1-F(x) ] (92) 
where f(x) is identified with 
ee i 2 
f(x) = &« exp(-x’) | exp(t ™)dt (93) 
O 


because this function has the same series expansion. Thus, 
finally, 


ei. 2 oe : 2 ae —, 
G(x, )=- O Sete if Xo exp(-x_,) J s=1,e xX, = y 


(94) 
What has been presented here is the linear theory for 
scattering. Appendix B gives calculated example curves of 
spectral distributions of plasmas with drift velocity W 
(W is smaller than Whee? the Grailtic abeveloci ty Jat whichmnon— 


linear drift instabilities are excited) W is in multiples 


L 
of the electron thermal velocity (3)? Ves a& is Oe the 
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Ee 


parameter identified in equation (30) and ¢€ = a - Each 
qe Be 
Ome 10° 


between W and k with ie = lwieos¢. The curves with the most 


eee 


spectrum contains 11 curves? from angles @ = 0, 


; T 
pronounced side maxima correspond to @ = ane the only symme- 


T ; : 
tric curve corresponding to @ = 2 with the horizontal tan- 


gent at y = 0. y is of course the dimensionless frequency 





, and extends over ly | < 3.5 representing the ion line. 
vi 
Relative maximum intensity I 1s given for each, since inten- 


sities vary strongly. Intensities are normalized so that 
| \a 12) + | iN 2 
1-G, (y)|" exp(-x') + 1G (x')I” 1 exp(-y*) 


TS ea ee er a a (95) 
| 1-G_ (x') = G.(y) | 





5) 
From the report by Theimer; by permission of 


Professor F. Schwirzke. 
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rit. SaEeARATUS 


The experimental setup itself requires a myriad of 
equipment. The schematic of figure 1 shows the equipment 
involved and its layout. The steady-state plasma facility 
at the Naval Postgraduate School is ideal for a systematic 
investigation of the laser light scattering plasma diagnos- 
tic technique. The parameters of plasma density and elec- 
tron temperature may be varied over wide range in the sys- 
tem. The plasma, whether argon, nitrogen, or helium, is 
produced through a hollow tantalum cathode and is acceler- 
ated to a disc anode between voltages of around 40 volts. 
The beam is columnated at the cathode by an arrangement of 
wound electromagnets producing fields of several hundred 
gauss. The plasma density may be controlled by varying the 
gas flow at a needle valve, altering the background pres- 
sure with a 10 inch diffusion pump, and by adjusting the 
field strengths of the columnating électromagnets. Plasma, 
and hence electron, densities; asy seen in the chee have 
a direct bearing on the scattered light intensities. Pre- 
paration of the system required the construction and place- 
ment of perpendicular window ports through which the inci- 
dent light entered and ‘the scattered light emerged. The. 
ports were approximately a half inch in diameter and were 
mounted on removable flames, as shown in figure 6. The 
iplane descrabed byethewpcerpendiewmlar incident mand scattened 


light wave vectors was arranged in such a way as to intersect 
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the plasma jet at a point within jaghalt inehnwor the tanta— 
lum cathode outlet, for two reasons. One, this is presum- 
ably where the beam and hence the electron density are ene 
greatest; secondly, little is known about the plasma para- 
meters at the cathode... the scattering technique should 
yield much information concerning these parameters. Further- 
more, laboratory space limitations and apparatus configura- 
tion demanded location at the spot chosen. 

The incident laser radiation was provided by the model 
130 high peak power Q-switched laser, manufactured by Optics 
Technology, Inc. The instrument consisted of a power 
supply/control console and a separate optical head, which 
incorporated a Fabry-Perot type interferometer and laser 
cavity mounted on a granite base. The Q-switch is a device 
which enables the laser to deliver extremely high peak power 
pulses varying from one to a hundred megawatts. As here em- 
ployed, it was a bleachable liquid system contained in a 
small cell with optically flat faces, which only needed to 
be inserted between the laser ruby and the interferometer 
for" operationvemFfurther funceticn*and™usé of "this" devaee 
will be discussed later. The laser crystal is a production 
quality, 60° orientation, three-eighth by three inch-.ruby 
crystal, emitting radiation at 6943 2. The crystal is op- 
tically pumped by two linear xenon arc lamps with three 
inch arc spacing and overall diameter of eleven mm, mounted 


opposite the crystal holders in a polished aluminum cavity. 
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Cooling is accomplished through a forced draft, water-cool- 
ed air system. The rear element of the Fabry-Perot inter- 
ferometer is a TIR prism, and a multi-element glass re- 
flector mounted on a micrometer driven assembly comprises 
the front element. Alignment of the interferometer optics 
is extremely critical for laSing action to occur. Align- 
ment may be accomplished in the following way. With the 

top of the laser cavity and the ruby crystal removed, care- 
fully position a cw laser so that its beam is in line with 
the axis of the ruby crystal holder pins and passing through 


the center of the exit hole as shown. 


CW LASER 


RUBY LASER | 





adjustable 
mount 


The task of aligning the cw laser may be more easily accom- 
plished accurately by placing a pen or a pencil on the ruby 
pins, holding it vertically, and adjusting the cw laser 
until the beam strikes the pencil in the center when placed 
on either pin. There will now be a multitude of reflec- 
tions on the face of the cw laser, one from the TIR prism 
and one from each of the glass reflectors in the microme- 
ter driven interferometer. Determine which one of these 


reflections comes from the ITR prism by placing a sheet of 
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paper directly in front of the prism; Once finding the 
correct reflection, remove the paper and hand adjust the 
prism until its reflection is as close to the cw exit aper- 
ture as Boeeipie. Then tighten its screws at its base to 
secure it. Next adjust the micrometer screws on the inter- 
ferometer until the multiple reflections made by its ele- 
ments merge on or near the cw exit aperture and prism re- 
flection. This accomplished, fringe patterns should be 
observed, and the ruby laser is aligned. 

The bleachable dye Q-switch device previously mention- 
ed permits the formation of giant pulses. The technique, 
usually referred to as passive Q-switching employs a small 
cell containing a solution of one of several metallo-organ- 
ic compounds known as phthalocyanines. In this case, a 
se-Lwtton—of Pace ain in methanomate used. THe. 61h- 
tion strongly absorbs light of the ruby Evaeeeney, and this 
absorption prevents net amplification Sfelignteunt: isa 
greater proportion of Gr?" gens have been pumped into the 
upper energy state than is normal. The energy input con- 
tinues until the eos tigers in the ruby overcomes losses 
due to absorption in the cell. To this point the laser 
smits coherent light weakly. A small amount of the light 
bleaches the cryptocyanine which then becomes nearly per- 
fectly transparent to the ruby light. . At the same instant 
there is a large net amplification and a giant pulse con- 


taining all the stored ruby energy develops rapidly. 
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Following the pulse, the solution returns quickly to its 


absorbing state. Below is a typical phthalocyanine absorp- 


tion curve: 


absorption 


9000 9500 6000 6500 7000 7200 7500 
wavelength 


Bleaching of the metal cryptocyanines is effected by the 
Saturation of the absroption at the ruby frequency. That 
is, the solution is in the bleached condition whrn the pop- 
ulation densities in the upper and lower states of the 
phthalocyanine molecule (S',S in the diagram below) are 
roughly equal. Because of the large cross-sections for 
capturing ruby light only small numbers of cryptocyanine 
molecules are involved. Consequently saturation occurs 


rapidly and relatively few quanta from the laser beam are 


required. 


Al 


excited triplet T’ 


first excited 
singlet state 
Ss’ 






triplet state T 







ENERGY LEVEL 
DIAGRAM FOR 
PHTHALOCYANINE 


14400 cm*! 


ground state 


The solution itself is prepared from the concentrate 
by diluting one part cryptocyanine to ten parts methanol. 
The cell is filled almost to capacity with methanol and 
nine drops of the 1:10 solution is added to this. This 
appears to give optimum Q-switching without requiring too 
high a threshold. 

Calibration of the laser output, that is, determining 
output energies and pulse widths (and hence output power) 
involves the use of a PIN photodiode (See Appendix C) and 
a laser radiometer. The HPA 4205 is a silicon planar PIN 
photodiode with a response time of less than one nanosecond 

2 


and a frequency response from dc to 10° hz. The diode was 


wired according to the following diagram. 
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PHOTODIODE 


20 = — 


“a 


The photodiode is useful in conjunction with a fast-rise 
oscilloscope, such as a tektronix 519, for recording laser 
pulse shapes. 

Laser output energy was measured using a Westinghouse 
model RN-1l Radiometer which is, more specifically a calor- 
imetric bolometer. It is designed to trap the energy beam 
from the laser in a bundle of fine, insulated copper wire 
which functions both as a calorimetric m ass and as a bolo- 
metric element. The change in resistance of the wire ele- 
ment 1s proportional to the energy absorbed and is almost 
entirely independent of the distribution of the energy with- 
in the unit. The bolometer element is part of a convention- 
al Wheatstone bridge circuit and the change in balance of 
the bridge is a measure of beam energy. (See figure 4.) 
When used in conjunction with an electronic microvoltmeter, 
laser energy may be read directly. When a microvoltmeter 
reads the imbalance in microvolts, the reading is directly 


Peoportioned to the incideémielaser pulse energy, 


fll ane io) aaa Laser output in joules 


in microvolts 
| (94) 


43 


The schematic of the calibration setup is as follows: 





beamsplitter 


ruby laser —_ 


oscilloscope 


The following are photographs taken on a Tektronix 
545A oscilloscope set to trigger as the pulse was initiated. 
The time scale on the three non Q-switched photos was 504u 
s/div., and on the fourth, Q switched pulse 500 ns/div. or 
O.5us/div. This gives an idea of the effectiveness of in- 
creasing lasing power by greatly shortening the pulse 
width. The Q-switched pulse width appears to be of the or- 
der of 200ns. The average pulse length of the non Q-switch- 
ed pulse appears to be of the order of 250uUs. Since power 
is energy per unit time, this shows an increase in lasing 
power by a factor of 1250. More specifically, a typical 
low power pulse will put out approximately 0.5 joules. For 
a pulse width of 250uUs this comes out to two kilowatts. 

The Q-switched pulse width of 0.2usec. the power is now 
2.5 megawatts, a substantial increase. The first three 
photographs show 3 pulses each; below them is recorded 
the average energy output of the three for a given input 
energy from the laser power pack. This was the same for 


each trace on a given photograph. 


Ad 


Non Q-Switched Pulses 
Sweep Speed = 50 Us/div. 


Input 1100 joules 


aidan aS oh Se 


BRREREE ss = Output 0.19 joules 





feet 1225 joules 
Seeput 0.38 joules 
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Non-Q-Switched 
Sweep Speed 50 Us/div. 


Input 1300 joules 
Output 0.55 joules 





For these first 3 photos, a plot of input versus output 
energy is made, establishing a linear relationship above 
laser threshold. 
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The fourth photo is included as a fine example of a 


Q-switched laser pulse. Sweeptime here is 0.5 wus/division. 





The scattered light is analyzed with the Baird atomic 
modified Eagle Mount three meter concave Grating Spectro- 
graph. The spectrograph is a Rowland circle type mounting 


for which the focusing conditions are: 


r = RcosA A = Angle of incidence 
r’= RcosB B = Angle of diffraction 
r = §lit to grating distance 
r' = Camera to grating 
distance 
R = Radius of curvature of 


the concave grating. 





N = Normal 
G = Grating 
S& = Sit 

C = Camera 


AT 


The Rowland circle has a diameter equal to the radius 
of curvature of the grating blank. The grating is mounted 
tangent to the Rowland circle. 


The usual grating equation applies. That is 


nA = d(sin A + sin B) (95) 
n = Order number 
X = Wavelength in Angstroms 
d = Grating spacing in Angstroms 


A and B are of opposite sign if they are on different 
Sides of the normal. Spectrograph controls include toggle 
switches for the grating, the plateholder (on which the 
photomultiplier tubes are mounted) and focus. 

Rotating the grating changes the angle of incidence 
A (and hence B). ‘The Veeder-Root counter is geared so as 
to indicate the angle of incidence, subtract from 1000. 
The focus toggle increases or decreases the plateholder to 
grating distance, the tilt control moves the plate (mount- 
ing the photomultipliers) in and out.about a vertical axis 
through the center of the slit. The latter two controls 
are not as important when using the photomultiplier as op- 
posed to a photographic plate when a careful focus is 
required. 

There existed no calibration chart of the instrument 
to predetermine the position of the photomultipliers for a 
specific wavelength. This must be accomplished experiment- 


ally, using a calibration source of Known wavelength, for 


Ag 


example a helium lamp. Helium has two lines close to the 
6943 A ruby wavelength, one at 6730 A and another at 7108A. 
Sighting visually through the plate aperture, these may be 
readily identified. With a slit cut in dark cardboard, the 
position at which these lines strike the aperture may be 
marked. lLinearly interpolating between the two wavelengths 
gives a rough position for placement of the photo multi- 
plier. A fine adjustment may be made when the whole system 
is aligned. The alignment procedure will be described 
shortly. For purposes of this experiment the proper grat- 
ing setting was found to be 991.2 (that is, with the angle 
of incidence A equal to 8.8 degrees.). The dispersion may 
also be calculated with some ease. This will become of 
great use when actual spectral distributions are analyzed 
carefully. This relation may be derived from the basic 


grating equation. 
nrX = d(sin A + sin B) (95) 


r = £(sin A + sin BP) 


letting d be a constant 


s+ ) = = cos B (96) 
a i 


or with A constant 


Ar 
dB 
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df = rdB (97) 
where r is the focal distance and 
is the increment along circumference 


for dB increment. 


QB _ L 
ak r 
(98) 
di _ gd | @B 
adi dB ak 
dk _ a ea 
re ie (99) 
bug 
r = Rocos B (100) 
a” mw cose D =eactual gratingespacimgs 
Dee dies Te) cose 
d = effective grating spacing. 
2. 
df mR on) 
d = Docos B (102) 
therefore 
an . D cos B 
qh nR EOE? 
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Numerical calculations for dispersion may now be made, let- 
ting D be the spacing for 15000 groove/inch grating and 
R = 3 meters at the 6943A~ wavelength and an angle of inci- 


dence A (determined experimentally) of gus, 


80 
i a 254mm LOA O 
De — T5000” inches x am x aoe = 16933A 








for the first order spectra 








X = d(sin A + sin B) A = 8.8° 
Sin B ae - sinA 
4 1 

_ 6943 Oo 

= aig053  “Winaeee 

= ,41003 GAN al oZ > 

= 2 ee 

cos B = .96638 
dX _ Dcos B _ (16933) (.96638) 
AL Re bane 10° 
O 
= 20280 = 5,458 A/a 
3 x LO 
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The photomultiplier tubes were RCA 7102 tubes with 


Sensitivity and response curves given below. 
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9000 7000 9000 11000 
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Supply volts 


The tubes were operated from a number of power sup- 
plies, the best of which appeared to be the Power Designs 
Model HV 1565 regulated dc power supply. Optimum signal 
to boise ratio was obtained at a voltage setting between 
1450-1490 volts (actual voltage as read on the voltmeter 


was closer to 1350 volts). The output signal was read 
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across about a five megohm impedance for Qtearer signal 
strength. This was of course at some sacrifice to response 
time, as will be seen later. 

Before discussing eHeiamoneane alignment procedure, 
the optics used to transfer the scattered light from the 
plasma to the spectrograph will be aese@ined: The arrange- 
ment was a periscope type as shown in figures 1 and 5. ‘Two 
converging lenses were employed, one focusing at the plasma, 
the other refocusing the plasma light at the exit of the 
periscope. A quartz condensing lens on the optical bench 
of the spectrograph refocused the light on the slit. An 
arrangement of rigid cardboard packing tubes and aluminum 
mirror holders formed the remainder of the device. The po- 
Sition of the periscope exit (on the plane of which the 
plasma beam is focused) is determined by using the grating 
lamps and quartz condensing lens in the following way. The 
grating lamps are hung over the front face of the grating 
holder and switched on from a control on the external power 
supply. The quartz condensing lens is placed in its holder 
On the optical bench and left free to move along the length 
of the bench. The image of the grating lamps may be observ- 
ed by holding a white sheet of paper behind the lens. Mov- 
ing the paper away from the lens will eventually focus the 
grating lamps in two distinct images. Moving the paper 


further away, the images will merge, and it is at that 
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point that the source, in this case the exit of the peri- 
scope, is to be placed. The lense is now secured and the 
periscope is mounted in its proper position. 

Because of the distance it was necessary to convey 
the lights and because of the small image of the lenses, 
alignment was extremely critical and somewhat difficult to 
obtain. Accurate laser-plasma-periscope-spectrograph align- 
ment may be obtained while the plasma system is at atmos- 
pheric pressure and the tantalum cathodes are being replaced. 
With the cathode flange removed, remove also the left win- 
dow flange. Place a wooden slate across two stationary 
bolts of the cathode flange in such a way as to have the 
slat horizontal and just below the center of the vacuum 
chamber (and hence below the plasma beam axis). The care- 
fully insert a rigid aluminum or stainless steel rod down 
the chamber, supporting one end inside the annular disk 
anode and the other end on the wooden slat. The slat 
should be of such thickness that the axis of the rod it 
Supports is parallel to and just below the chamber axis. 
This arrangement will now support the mirror proble to be 
inserted from the left window opening; the probe will act 
as an artificial "Scatterer” so that proper alignment both 
Or the Optical transfer system and of the grating’ may Dbe.ob- 
tained. The mirror probe consists of a straight wooden 


dowel with one end cut to 45°. To the 45° cut is applied 
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a small cutout piece of polished stainless steel to act as 
a reflector. This may be glued on and a mark may then be 
made on the opposite end of the dowel to indicate the pro- 
per vertical mirror position. The dowel is then inserted, 
one end supported on the rod arrangement described above 
and the other end on a notched cardboard edge, taped to 
the window opening. The supporting rod may vary the position 
of the reflector along the chamber axis without significantly 
changing its position perpendicular to the axis; pushing 
or pulling the mirror rod itself, varies its position per- 
pendicular to the chamber axis. 

The next step is to place a strong light source such 
as a mercury lamp before.the grating (inside the specto- 
graph) and aimed back toward the slit. The Mercury beam may- 
be focused at the slit using its own lens for more intensi- 
tys eWith theebaboratery dawkened; carefulidyeadguciteethempen- 
iscope arrangement until the light from the Mercury falls 
on the supporting rod, which has been located parallel to 
the vacuum chamber center axis. This done, the mirror prop 
may be adjusted so that the light falls on its reflector. 
From the end of the mirror proble, and peering into the 
chamber through the left window mount the probe may be ad- 
justed carefully so that its light falls on the half inch 
laser window. The axis of the mirror should be vertical, 
that is, reflecting straight up. With the Mercury light, 
(which may be easily discerned) shining on the half inch 


laser entrance window, the system, save the laser, should 


Sie: 


now be in perfect alignment. The next problem is to align 

the laser so that its artificially "scattered" light may be 
seen by the grating; and when this is done the grating must 
be rotated until the laser frequency is centered on photo- 

multiplier tubes at the spectrograph exit. 

The laser is mounted on an aluminum vertically adjust- 
able stand, the laser itself containing tilt adjustment on 
its base. Alignment along the horizontal plane may be made 
more easily with a cw laser by raising the axis of the cw 
laser, the pulse laser may replace along the same axis and 
trial and error adjustment made until alignment in the ver- 
tical plane is achieved at will. Once aligned the laser 
pulse should be seen in the slit when fired. 

The grating may then be rotated to the proper wave- 
length. The phototubes and grating have already been rough- 
ly adjusted according to the helium line procedure described 
above. First the PIN photodiode is positioned with an arran- 
gement of beamsplitters (microscope slides) to analyze the 
laser pulse. The signal from the PIN photodiode is then fed 
to the upper trace of a dual beam oscilloscope, and the sig- 
mae £rom the photomultiplier to the lower trace. Upper 
trace gain is lowered because of the intensity of the laser 
pulse, but the oscilloscope arrangement is as shown in 
figure 14. The grating is rotated by 0.1 degree at a time 
in order to find the optimum position. Below are the traces 


illustrating how this was accomplished. The upper trace is 
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the pm tube signal and the lower trace is the laser pulse. 
Beside each is the respective grating setting and its cor- 


responding angle of incidence. 


991 28. o 
991.1 "eee ) 


991/50 cme) 
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991.5 (8.8°) 


oo. gee. 7) 


991.2 (8.8°%) 





Sweepspeed was 50 us/div. Both traces were triggered by 
the laser pulse, the lower trace of each group. The laser 
trace gain was 0.5 V/div. and the pm tube trace 25 V/div. 
All were non Q-switched pulses. Note the difference in 


pulse shapes. As was mentioned before, this was due to 


Seo, 


the large resistor in the signal line (about five megohms) 
which considerably lengthened the rise time of the tubes. 
This of course, 1S no consequence since relative intensi- 
ties are what is important here and in analyzing the spec- 
tral distribution of the scattered light. 

The system, now constructed and aligned, may be em- 
ployed in the attempt to discern a scattered light signal 
and analyze the spectral distribution of the scattered 
light. The few attempts thus far have been unsuccessful. 
Some combinations of detecting tubes and oscilloscope am- 
plifiers have been tried to increase the signal to noise 
ratio. Photographs have been taken, an example of which 


is below, 
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The photograph on the bottom was triggered at the 
scope with no laser pulse in order to provide a noise refer- 
ence from the pm tube. Tube voltage was set at 1490 (Ac- 
tually 1390 on the voltmeter). Sweep speed was 50 us/div. 
and the background pressure in the vacuum chamber was 
8 xX io mm Hg. Oscilloscope gain here was 0.05 volts/div., 
as shown in figure 14. If a signal is present, as a few of 
tae traces such as the elle doove scem tO indicate, it is 
virutally impossible to distinguish it from the noise. 
Possible improvements will be discussed in the conclusion. 
Attempts at focusing the laser beam seem to have no 
effect, but this is almost certainly due to the poor quali- 


ty of lenses available. 
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Iv. CONCLUSIONS 


The scattering apparatus has been erected and aligned. 
The attempts which have thus rar been made, to discern the 
scattered light and analyze its spectrum have been unsuccess- 
ful. The plasma parameters specifically the density, affect 
scattering intensities as much as many varied instrument 
parameters affect the experimenter's ability to detect the 
light. All oscilloscopes, save a very few are in dire need 
of realignment and recalibration, as are the plug-in units 
that accompany them. When the scattered light is finally 
detected up to seven photomultipliers will be needed to si- 
multaneously analyze the spectral distribution of a single 
pulse; knowing the dispersion of the grating, the photo- 
multipliers, connected to the spectrograph plate with short 
lengths pipe, may be accurately spaced around the central 
frequency to observe the whole spectral distribution of the 
scattered pulse. These should be of the RCA 7102 type. 
The 7102 tubes now available appear to have lower signal 
to noise ratios at sensitivities required that one would 
expect. Improvement of that ratio could be of significant 
importance in detecting the scattered light. Increasing the 
plasma density by purging more gas into the plasma system, 
increasing the field on the cathode and anode magnets and 
by controlling the main diffusion pump may well be the 


answer sought. Furthermore, the laser system available, 
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while powerful, may not be quite what is needed, though 
focusing of its beam through a lens of better quality may 
help. A more energetic system is now on order. It is 
interesting to note previous experimental attempts to study 
at light scattering and the problems mentioned in these 
efforts. One group, Thompson and Fiocco at MIT measure the 
light of a laser scattered from an electron beam. The elec- 


tron density was about 5 x lige em, and the ratio of 


scattered to incident intensity was of the order of ow 
Applying their technigue to a low electron density plasma, 


n Ot? cm>, they improved the ratio of scattered to in- 


cident power to Meme 
The best scattered to incident ratio came from E. 
Funfer, B. Kronast, and H. J. Kunze, who measured light 


scattered from a high density pinch plasma (n Mio). 


Their ratio came to Tre 

Distinguishing the scattered light from background 
plasma radiation is a considerable problem. One possible 
solution that may be tried is to eliminate the plasma ra- 
diation with polarizers, since at 90° the scattered iLiel ake 
will be polarized, whereas the plasma background radiation 
will not. One other wasy of eliminating plasma radiation 
would be with the use of a differential amplifier which will 
differentiate the scattered signal from a pair of pm tubes 
one seeing the scattered light and radiation, the other only 


radiation. This technique was attempted, but the character- 


istics of the two pm tubes available were not close enough 
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for comparison. Furthermore, the differential amplifier 
plug-in units on the oscilloscopes were not functioning 
properly. 

The above has summarized experimental progress to 
date. This progress should enable future experimenters to 
resume development of the experimental technique of laser 


light scattering as a plasma diagnostic. 
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APPENDIX A. DERIVATION OF THE BOLTZMANN EQUATION® 


Each electron is represented by a point in six dimen- 
Sional electron phase space. The coordinates of the phase 
space are xX, Y, 2, U,, Wags u,. Because each electron has 
a velocity, the representative point for each electron must 
change its position coordinates with time. In addition, 
any change in the velocity of the electron will also cause 
the representative point to move in electron phase space. 
The velocity of an electron will change if a force acts on 
the electron or if the electron collides with any other 
particle. 

Consider a fixed differential volume in electron phase 
Space dV = dx dy dz du dy du, (the "electron" subscript 
will be omitted). An expression will be written for the 
number of representative points in this volume as a func- 
mromwOr ©, Uy and. ts egethat 1s;,e 1 willl besa sicontinua ty” 
equation for the representative points in electron phase 
space. The increase with time in the number of electron 

| 8(f£_dv) | 
mpoints" in AV, —y, 1s caused by three factors: 

1. The drift of points in position coordinates due 

to inherent velocity. 


2. The drift of points in velocity coordinates due 


to electric and magnetic forces. 


6 
From Uman, p.34-36. 
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3. The spontaneous velocity change without position 

coordinate change due to collisions. 
When these three terms are evaluated, the Boltzmann equa- 
tion for electrons may be written. Note that the effects 
of external gravitational forces and inteparticle nuclear 
forces are omitted. 

First determined is the increase with time of the 
number of representative points in dV due to drift in po- 
Sition coordinate space. The number of representative 
points which cross the surfaces of dV in time de because 
of inherent velocity must be calculated. The velocity of 
apeointeinethessedirection at x, YZ, u; Me u, is uy. 
The number of representative points in a volume of "base 
area" dy dz du, oh du, and height u,dt will enter dV in 
time de through a face normal to the x axis at x, y, Z, 


u fered, use Thsenunber@is 
3 (don) 
| | dt =f Gea eotayvdzdusaumdau (1) 
dt ~ eee x x Yy 2 


The change in the number of representative points in dV in 
time dt due to the exit of points through the equivalent 


surface at x+dx is given by 


Ol(cemii.) 


— 
S(d'N) | — = , _. eae. | 
[ ae er: fue rye cia) clit OY Cae 


(2) 
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Therefore the net rate of change of points in dV due to 
velocity in the x direction is 


Swi y 


6 6 
a | o(d N = 
[ Ot +{ “Gu ou Ox a (3) 


x+dx 


Thus the total rate of change of points in dV due to in- 


herent velocity may be written as 


6 
| 2{gw) ct ie = - (aay : Vf.) dV (4) 


Now consider the drift of representative points in velocity 
coordinates due to electric and magnetic forces. The 
"velocity" of a oS ee in the u,, direction in 
velocity coordinates is = x) = C sx) Using this 
"velocity" the net influx of points through the two surfaces 


perpendicular to u, may be calculated: 


6 - 
| ey + E34 7 ~ Ley © a, )] av 


(5) 
The extension of this expression to include the flow through 


the surfaces perpendicular to u_ and u, follow 


(aon) 


agin] --(B-™ 8) a : 
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At this point, the time rate of increase of represen- 
tative points in dV is equated to the influx from the two 
drifts (neglecting coliisions) and the Boltzmann equation 
is obtained: 

1 


Oe 
fe > dv = - € © Vr ) dv - mo °Vut )av 
t —e = m = 


or (7) 
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APPENDIX C 
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a Me a 





ULTRAFAST LOW NOISE 
SILICON PIN PHOTODIODE 


The hpa 4205 is a silicon planar PIN photodiode for 
etection of light in the visible and near infrared regions. 
esponse to blue and violet light is unusually good for a 
ery low, dark current silicon photodiode. 

Speed of response of this diode is less than 1 nanosecond. 
aser pulses shorter than 0.1 nanosecond may be observed. 
he frequency response extends from DC to 1 ge. 

The low dark current of less than 150 picoamperes en- 
les detection of very low light levels. The quantum detec- 
on efficiency is constant over six decades of light intensity, 
roving excellent dynamic range. 

An extremely fast photon coupled pair can be obtained 
y using the hpa 4205 in conjunction with the hpa 4107 
nfrared light source in an identical package), for con- 
mnient use in card and tape readers, encoders, and similar 
»plications. 


OPTICAL CHARACTERISTICS 
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ELECTRICAL CHARACTERISTICS (25°C) 
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TECHNICAL DATA 15 FEB 66 





Test Conditions 

0.7 Sselectrons,phatomeesee.aeeeee:...... Mae ., Wesenter = 20) V 
0.25,.4,/mW /CM?*... 5a. mee 8 .....cesutl:. ee R, = <1Mon 
0.125 pA/pWw 
0.5 x 1073 cm? 
0.010” diameter 
Less than 1 nanosecond 

Test Conditions 
Fabs) 0101 0% appa 7 Vr = — 10 volts 

plover Vacsci (holy et (ele: - coe... | aaa ee Vp = 1 volt 

ace eeuwcemsnne ssa eesiceciescca alc ea ceecemececneeaacecnencs sce eceeccaceecoke Va a 10 V 


See OPP OOS OHH SO SSO HK SeeeerFFeFFeeteesFeseses Feet 2Geen ==2 528 


MAXIMUM RATINGS 
Steady Reverse Voltage 
Peak Inverse Voltage 
Power Dissipation 


50 volts 
200 volts (see Note) 
5SomW T= 25°C 


qP2 


Note: Exceeding the PIV will cause permanent damage to 
the diode. Forward current is harmless to the diode, within 
the power dissipation limit. For optimum performance, the 
diode should be reverse biased at between 5 and 20 volts. 
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Figure 1. Spectral Response 
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Figure 6. Relative Directional Sensitivity of the PIN Photo- 
diodes. 
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Figure 8. Typical Output Characteristics. 


Cui 


OR TTKE 
PAR 


RELATIVE RESPONSE 





WAVELENGTH (MICRONS) 
Figure 9. Spectral Response. 


(pf) 


CAPACITANCE 





TEMPERATURE (°C) REVERSE atk . | 
Figure 10. Typical Capacitance Variation with Applied 


Figure 7. Dark Current at — 10 V Bias vs. Temperature. Voltage. 





APPENDIX D 


9d09s 


19S 0 


INIWIY dX] 





SNIMILIVIS YISV) 


aN 
apoipoyoyd® 
\ 
\ 


JY!) pasaqjeors 


at SU 9 jw 


adoasiiaed 


? 


| angi) 


S139} jdsweag 


82 


IOOHLVI 


dwnd woisnys!p,.01 


t 


}aude a ib —” - 2 ae 


ee Pacer 4 
YP ew SmPtnwaerarreEaasaaats, 
~~ @ Cee Se aeeenrnwae ear 2a t= 
i wa oe 2 Pee» iT 12 easwvarrat 
P—— | | ome 





| Ya on om 
= ry ry. 
WHAUUY 
a yno oul 


13]eM 
dui|009 


tj 
ee 


| 





omits ~ sedis tee Be. fy 


(ES!) es, 


Z aindi} 


JOONV 


83 





REMOVE TRIGGERING 


LASER HEAD 


84 


head 


contro! unit 


25 


on—off switch 





on 


os if} = 


6V 


CIRCUIT SCHEMATIC OF RN-1 RADIOMETER 


fig.4 


85 


‘ili 

















ee ee 
e2cer2eace® 


- oom 
rpe@c@ercersarevaree 


ore ae? 


> — 
o . ed oo . 










a ee 





ay 


Po 


§ 
?- 


atin 


as ay 


87 





88 





Sig, 


S1ayijdswmeag 


—- 





90 


ewseld 


apoue 
ISI 
lejnuue 





grayscale 
7% 


e ¥ 
ee, ? —* 


Siouuod wandea pue aiyewayos saquieyg,... |<) 


———— pore a a 


—te 


—_ 





sual 


dUISUapU0d 





a3 








oscilloscope 


‘s C settings 


HORiZ POSITION 
HORIZ DISPLAY 

4 mt B eu;, 

a: xe x2 











see 
N10 *s 


yt 


ua 


x 

































TIME BASE A » -_ 
TYPE 2) TIME- BASE “UN IT 7” — - - 
i 
TWeErE —— TYPE 22 TIME-BASE UNIT 
SLOPE COUPLING wrote Oman 
. etem etiam 
sg. 5 time Ont an 
Bicuetan? ts 
VARIABLE “Galt vane 1 Gate 
TIME/CM san etoni TIME CM 
ee ee 


FUNCTION 





StASHITT Sswerrs One 


toe 14cm 4 
BNO tee | 
C8 bt At - 








j SAWTOOTH SINGLD SwEe a 
¢ SEC Pula TO 
eases 
ear 
Cause arty 
Com raOmek ine Teetsonig sme 
SERIAL PORMERO CRIGONM U5 & ; SERIAL COOMAND COIGEM USA 





ies 





UPPER BERN tsi uta vette amr : ave LOWER BEAM 


C emeneen eeeemiina teed 









TYPE H CALTN-VSN- Trawms-y F TYPE Z PLUG-IN UNIT 
. WIDE-BAND $ eee tinve. scat atmo 
PLUG-IN UNIT | _woreane | PLA MA Fa Smet eS 


3 aav. 06 v/Cm OC COVeLED SERTICAL 
POLAMP QISETIME «= 3 we SEL 


FOSITION 






e@ 


INPUT A 









GAatw ADJ. 
pey e 
INPUT 
SELECTOR 
Oc my, 
y Diff ¢ A-B 
- . OiFe 
VARIAGLE 
VOLTS/CM 


GAIN VAR ATTEN 
acsust BALANCE 


©, AS 
By 


we Onteoe OS 4 


COak SE 





VOLTS/CM 


UFe Pet tot 


e 4 


wmtasiseale® CaLieeates 





SERIAL 





~ 












PORTLAND, OREGON, U.5. A. Heeseoreia tee 





PORTLAND, OREGON, U.S.A. 


oS 


ie 


2.3 


L@e 


der 


3. 


14, 


TS. 


ts 


BIBLIOGRAPHY 


Baird-Atomic, Inc. Manual, The Three Meter Grating 
Specteogmapn, EG3. 


Chan, P.W. and Nodwell, R.A., "Collective Scattering 


of Laser Light by a Plasma", Physical Review Letters 
16@ED. 12278 1966. 


Fiocco, G. and Thompson, E., "Thomson Scattering of 
Optical Radiation from an Electron Beam", Physical 
Review Letters, 10, p.89, 1963. 


General Dynamics Corp., General Atomic Division Report 


1PP 1/48, Light Scattering and Drift Instability ina 
Plasma by Otto Theimer, translated by Fred Schwirzke, 


p.2ff, July, 1966. 


Greim, H.R., Plasma Spectroscopy, pp.180-184, McGraw- 
Hime Beok Com, 1904. 


HPA Corp. Manual, The PIN Photodiode, 15 Feb., 1966. 
Huddelstone, R.H. and Stanley, L.L., (Editors), Plasma 
Diagnostic Techniques, pp. 594-596, Academic Press, 
Them, 1968. 

Institut Fur Plasmaphysik, Report IPP/6/9, The Laser as 


a Tool for Plasma DLagncs tres. by W.H. Kegel, pp.2ff, 
Oct., uae . 


Optics Technology, Inc. Manual, The Model 130 High Peak 
Power Q-switched Laser, Apr., 1966. 


Salpeter, E.E., “Density Fluctuation in a Non-equilib- 
rium Plasma", J. Geophys.Res. 68, p.4870, 1963. 


Salpeter, E.E., "Electron Density Fluctuations in a 
Plasma", Physical Review, 120, p.1528, 1960. 


Schwirzke, F., Personal Communication. 


Smith, W.V. and@iiserokin, ©.Pyemeeeaser, pp.lsi-1>2 
McGraw-Hill Book Co., 1966. 


Uman, M.A., Entroduction sto Plasma, ehysiae, po, 3336, 
McGraw-Hill Book Co., 1964. 


Westinghouse Electric Corp., Instruction Manual, The 
Model RN-1l Laser Radiometer. 


96 


INITIAL DISTRIBUTION LIST 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0212 
Naval Postgraduate School 
Monterey, California 93940 


Naval Air Command (Code AIR-310) 
Department of the Navy 
Washington, D.C. 


Assoc. Professor Fred Schwirzke (Code 61Sw) 
Department of Physics 

Naval Postgraduate School 

Monterey, California 93940 


Professor N. EL. Oleson (Code 610e) 
Department of Physics 

Naval Postgraduate School 
Monterey, California 93940 


LTJIG Alfred A. Pease 
a2) M2) 33rd Road 
Arlington, Virginia 22207 


No. Copies 


20 











Unclassified 
Security Classification 


DOCUMENT CONTROL DATA-R&D 


(Security ciassification of titie, body of abstract and indexing annotation must be entered when the overail report is classified 
2a. REPORT SECURITY CLASSIFICATION 


Unclassified 


2b. GROUP 


Laser Light Scattering As a Plasma Diagnostic Technique 


1. ORIGINATING ACTIVITY (Corporate author) 






Naval Postgraduate School 
Monterey, California 93940 











3. REPORT TITLE 










4. OESCRIPTIVE NOTES (Type of report and inclusive dates) 


Master's Thesis; ne 1969 
5. AUTHOR(S) (First name, middle initial, last name) 


Alfred A. PEASE 


6. REPORT DATE 7a. TOTAL NO. OF PAGES 7b. NO. OF REFS 
| eaupe 1969 98 16 


6a. CONTRACT OR GRANT NO. 98. ORIGINATOR’S REPORT NUMBER(S) 

















6b. PROJECT NO. 





9b. OTHER REPORT NO(S) (Any other numbers 
this report) 





that may be assigned 








10. OISTRIBUTION STATEMENT 





| his document hes «en approved for public 
eleas 3] stribution is unlimited. 









n 








11. SUPPLEMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 













Naval Postgraduate School 
Monterey, California 93940 





13. ABSTRACT 





The effort here was concerned with initiating a long term experi- 
mental study to develop the method of laser light scattering as a 
plasma diagnostic technique. This was done with a ruby laser and a 
grating spectometer coupled with sensitive detectors, at the steady 
state plasma facility at the Naval Postgraduate School. Although 
attempts to detect the extremely low levels of scattered light failed, 
excellent progress was made in determining existing apparatus limita- 
tions and in recognizing which system parameter may be adjusted in 
order to detect the scattered light; this will enable the future 
experimenter to progress with a minimum of overlap. 









FORM : 
DD 1 NOV 21473 Le Aoes |) Unclassified 
S/N 0101-807-6811 Security Classification 


=) 5, 


‘A-31408 


Unclassified 
~ Security Classification 


KEY WORDS - aie 


Laser Light Scattering 

Plasma Diagnostic Technique 

Scattered Spectral Distribution 
: i . 8 


DD Nore 14 7 (BACK) Unclassified 


S/N 0101-807-682}3 Security Classification A-31409 


100 


























es 
Laser light scattering as a plasma diagn 


